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Abstract Using axisymmetric MHD simulations we investigate how the overall 
jet formation is affected by a variation in the disk magnetic flux profile and/or the 
existence of a central stellar magnetosphere. Our simulations evolve from an initial, 
hydrostatic equilibrium state in a force-free magnetic field configuration. We find 
a unique relation between the collimation degree and the disk wind magnetization 
power law exponent. The collimation degree decreases for steeper disk magnetic 
field profiles. Highly collimated outflows resulting from a flat profile tend to be 
unsteady. We further consider a magnetic field superposed of a stellar dipole and 
a disk field in parallel or anti-parallel alignment. Both stellar and disk wind may 
evolve in a pair of outflows, however, a reasonably strong disk wind component is 
essential for jet collimation. Strong flares may lead to a sudden change in mass flux 
by a factor two. We hypothesize that such flares may eventually trigger jet knots. 



1 Jets as collimated MHD flows 



Astrophysical jets are launched by magnetohydrodynamic (MHD) processes in the 
close vicinity of the central object - an accretion disk surrounding a protostar or a 
compact object [1,2, 9, 20, 21, 24]. Numerical simulations of MHD jet formation 
are essential for our understanding of the physical processes involved. In general, 
simulations may be distinguished in those taking into account the evolution of the 
disk structure and others considering the disk surface as a fixed-in-time boundary 
condition for the jet. The first approach allows to directly investigate the mecha- 
nism lifting matter from the disk into the outflow [3, 10, 1 1, 14, 17, 18, 22, 24] This 
approach is computationally expensive and still somewhat limited by spatial and 
time resolution. In order to study the acceleration and collimation of a disk/stellar 
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Fig. 1 Example initial magnetic field distributions (poloidal magnetic field lines). Full and dashed 
lines indicate the direction of magnetic flux. Magnetic field parameters: Ifo.d = 0.01, —0.01, —0.1, 
resp. f ,* = 5.0, 5.0, 3.0 (from left to right). From [5]. 



wind it is essential to follow the dynamical evolution for i) very long time ii) on a 
sufficiently large grid with iii) appropriate resolution. For such a goal, the second 
approach is better suited [4-7, 12, 13, 15, 19, 25], allowing as well for parameter 
studies. The case of superposed stellar/disk magnetic field is rarely treated in simu- 
lations, still, the first model was discussed already in [23]. Simulations of a dipole 
with aligned vertical disk field are presented by [16, 18]. The stellar field has impor- 
tant impact on the jet formation process as enhancing the magnetic flux, adding a 
central pressure, and providing excess angular momentum for the launching region. 



2 Model setup 

We use the ZEUS-3D MHD code extended for physical magnetic resistivity (see 
description in [6]). The set of MHD equations considered is the following, 

$- + V-(pv)=0, VB = 0, — j-VxB, ?-Vx [ vxB-- rjj ) =0 (1) 

at c at \ c J 

+ V( jP + A4 )+pV<£- J -^- B - =0, (2) 

c 

with the usual notation [4-7, 19]. We do not solve the energy equation, but apply 
an internal energy e = p/(y— 1) of a polytropic gas (y = 5/3). Turbulent Alfvenic 
pressure p& allows for a "cool" corona. The turbulent magnetic diffusivity 77 (r,z',t) 
can be related to p& applying our toy model [6]. The 77 ~ 0.01 was chosen low 
and does not affect collimation. Diffusivity is, however, essential for reconnection 
processes. We distinguish setup DW (pure disk wind) and SDW (stellar wind plus 
disk wind) by choice of boundary and initial conditions. Model DW investigates 
different disk magnetic field and mass flux profiles [4]. Model SDW investigates the 
interrelation of the stellar magnetosphere with the surrounding disk jet [5]. 



du 

dt 
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Fig. 2 Collimation degree 
< £ > and power index of 
the disk wind magnetization 
profile Bars indicate 
simulations with time variable 
collimation degree (from [4]). 
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Boundary conditions: In setup DW we distinguish along the equatorial plane the 
gap region r < 1.0 and disk region r > 1.0. The magnetic field is fixed in time and 
is determined by the initial condition. We have chosen a power law, B p (r,0) ^r - ' 1 , 
and investigate different jj.. In setup SDW we further distinguish the star from r = 
— 0.5, and the gap from r = 0.5 — 1.0. Co-rotation radius and inner disk radius 
coincide. A Keplerian disk is the boundary condition for the mass inflow from the 
disk surface into the corona. Matter is "injected" from the disk (and the star) with 
low velocity Vi n j(r,0) = ViVK(r)Bp/Z?p and density p; n j(r,0) = r]ip(r,0). Typically, 
V; ~ 10~ 3 and rjj ~ 100 for stellar and disk wind, but could be chosen differently. 

Initial conditions: As initial state we prescribe a force-free magnetic field and 
a hydrostatic equilibrium p{r,z,t = 0) = (r 2 +z 2 )~ 3//4 . For model DW we calcu- 
late the initial field distribution from the disk magnetic field profile using our finite 
element code (see [4,8]). For model SDW the initial field is a superposed dipole 
plus disk field. For the disk component we apply the potential field of [6,19]. We 
prescribe the initial field by the magnetic flux distribution f(r,z) = /B p <iA, 

f (r,z) = fo,d - (Jr 2 + (z d +z) 2 - (z d +z)) + ^ /2 . (3) 

Certain field combinations are investigated, parameterized by the disk ffo <j and stel- 
lar magnetic flux ff) * (Fig. [TJ. 



3 Disk wind magnetization and jet collimation 

Simulations of setup DW were run for different disk magnetic field profiles B p (r,z = 
0) ~ and density profiles Pi n j(r,z — 0) ~ r~fp . In general, we find an increasing 
degree of collimation with decreasing slope of the disk magnetic field profile (see 
[4]). This seems to rule out launching models for collimated jets from a concen- 
trated magnetic flux such as e.g. the X-wind scenario. A steep density profile leads 
to a higher collimation degree, which is not surprising as the mass flux is more 
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Fig. 3 Poloidal magnetic field evolution during one example flare around r = 180. Solid and 
dashed lines indicate the direction of total magnetic flux of the supeiposed dipolar and disk mag- 
netic field components. Shown are time steps: 1760, 1790, 1810 (from left to right). 



concentrated just by definition of the boundary condition. A physically meaningful 
classification taking into account both density and magnetic field can be achieved by 
comparing the degree of collimation degree versus disk wind magnetization profile, 
a{r,z = 0) ~ Bl(r)r 4 p^v^{r)Q K (r) 2 , thus, a(r,z = 0) ~ r-Qp-Vp+W) = r ^ . 
The resulting diagram Fig. [2] shows a convincing correlation between the magneti- 
zation power law index jj, a and the average degree of collimation <£>■ The width 
of the (Ho - <£>)-correlation is due to further differences in the parameter space. 



4 Jet mass flux triggered by star-disk magnetospheric flares 

Simulations of setup SDW were run for aligned and anti-aligned orientation of 
dipole versus disk field and for different strength of both field contributions [5]. 

Independent of the alignment, the central dipole does not survive on the large 
scale. A two-component outflow emerges as stellar wind plus disk wind. For a rea- 
sonably strong disk magnetic flux a collimated jet emerges. If the overall outflow 
is dominated by a strong stellar outflow a low mass flux disk wind remains uncolli- 
mated. The best setup to launch a collimated jet from a star-disk magnetosphere is 
that of a relatively heavy disk wind and high disk magnetic flux. Stellar wind domi- 
nated simulations may give a high degree of collimation, however they collimate to 
too small radii. Stellar magnetic flux dominated outflows tend to stay un-collimated. 

In some simulations we observe reconnection flares, similar to coronal mass ejec- 
tions, typically expanding and reconnecting within 70 orbital periods of the inner 
disk. This is similar to [10], however, in their case reconnection is are triggered by 
time-variation of the accretion rate. In our case the reconnection/flares seem to be 
triggered by the evolution of the outer disk wind. Even for our very long time-scales 
the outer disk outflow is still dynamically evolving, thus changing the cross-jet force 
equilibrium and forcing the inner structure to adjust accordingly. The flare events are 
accompanied by a temporal change in outflow mass flux and momentum. Figure |4] 
shows the mass loss rate in axial direction integrated across the jet. We see two flares 
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Fig. 4 Axial mass flux in- 
tegrated along the upper 
z-boundary versus time. Note 
the change of mass flux of 
10 - 50% during the flare 
events. High mass fluxes for 
t < 500 indicate sweeping off 
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with a 10%-increase in the mass flux followed by a sudden decrease of mass flux by 
a factor of two. This behavior is also mirrored in the poloidal velocity profile. 

Considering the ejection of large-scale flares and the follow-up re-configuration 
of outflow dynamics, we hypothesize that the origin of jet knots is triggered by such 
flaring events. Our time-scale for flare generation is of 1000 rotational periods and 
longer than the typical dynamical time at the jet base, but similar to the observed 
knots. The flare itself for about 30-40 inner disk rotation times. 
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